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HEAT  TRANSFER  TO  A CYLINDER  FOR  THE  FREE 
MOIECUIg  FLCW  OF  A NON-UNIFORM  0 AS 


1.0  INTRODUCTION 

The  heat  transfer  to  a body  In  free  molecule  flow  has  been 
determined  theoretically  by  Tsien  (Ref.  l),  Stalder  and  Jukoff  (Ref.  2) 
and  others  (Refs.  3 and  4)«  Experiments  by  Stalder,  Creager  and  Goodwin 
(Refs.  3 and  5)  have  provided  some  confirmation  of  the  theoretical  de- 
velopments. 

These  investigations  have  been  confined  to  gas  flows  in  Max- 
wellian equilibrium.  The  present  report  is  a continuation  of  Reference 
6 in  which  the  aerodynamic  forces  on  a cylinder  were  determined  for  the 
case  of  a non-uniform  gas  (i.e.,  one  in  which  heat  flux  or  viscous  stresses 
are  present). 

It  is  assumed  that  the  right  circular  cylinder,  which  is  ori- 
ented perpendicular  to  the  gas  flow,  has  so  large  a thermal  conductivity 
that  it  has  a uniform  wall  temperature.  Possible  end  effects  are  neglected. 

The  results  depend  on  a knowledge  of  the  molecular  velocity  dis- 
tribution and  the  distribution  assumed  here  is  the  one  proposed  by  Grad 
(Ref.  7).  The  results  are  complete  for  the  monatomic  case  but  only 
approximate  for  the  polyatomic  case. 

For  the  sake  of  brevity,  much  of  the  general  discussion  given 
in  Reference  6 has  not  been  repeated  here. 

2.0  METHOD  OF  ANALYSIS 

The  general  method  of  analysis  is  the  same  as  ip  Reference  6 
and  the  energy  balance  for  a differential  area  is  written 


dEi  -*-dTL  ■=  dE a,  -hd' Ra,  + <AQ 


(2.1) 
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Assuming  no  radiation  and  using  the  themal  accommodation  coeffi- 
cient 

ot  = ol£l  - 

cLZ  i — dEw  (2.2) 


one  may  rewrite  the  energy  balance  for  eonveotion 

•i-  dQ  - dE-,  - dEw  (2.3) 


The  incident  energy  flux  per  unit  area  per  unit  time  may  be 
expressed 

00  *o  °o 

(n'z  v,z  + wa)  ix'f  d.u.'dv'dw' 


d£t 


nm 

2 


(2.4) 


- oc  > <3 


where  £ is  the  distribution  function  used  in  Reference  6, 
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Kxoept  for  the  molecular  velooity  components,  all  the  variables 
appearing  in  Sqs.  2.5  and  2.6,  are  treated  as  constants  throughout 
the  region  oooupiad  by  the  cylinder.  The  relations  between  velocity 
components  are  given  by 

UL  — — \jJ  cos  9 «f  y 'sinB 
v = — uJ  send  — v'cos# 
w * w' 
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The  energy  flux  per  unit  area  per  unit  time  due  to  diffusely 
reflected  molecules  is 
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and  since  d£^/dik  is  to  be  computed  on  the  assumption  that  all 
of  the  molecules  are  reflected  diffusely 
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Combining  Eqs.  2.3,  2.4*  2.8  and  2.9  and  integrating  over  the  cy- 
linder ATT  00  °°  010 
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In  a polyatomic  gas  there  are  internal  energies  of  rotation 
or  of  vibration  or  of  both.  To  calculate  the  incident  energy  flux 
in  this  case,  one  should  use  a distribution  function  which  includes 
the  distribution  of  rotational  and  vibrational  energies  (e.g.,  see 
Ref.  8).  If  the  gas  molecules  have 3 degrees  of  freedom  of  motion 
(thus, j - 0 for  a monatomic  gas)  the  number  of  incident  molecules 
may  be  written  in  the  form 


co  oo  c»  ao  06 


duu'  dv  'dw  ' d V t • • • dVj 


(2.11) 


- 4 " 


where,  for  the  sake  of  simplicity  in  illustrating  a point,  the  clas- 
sical mechanical  approach  has  been  adopted  and  a characteristic 
velocity  has  been  associated  with  each  degree  of  freedom  of  motion. 
Actually,  an  accurate  calculation  would  require  a quantum-mechani- 
cal formulation. 

Assuming  that  the  total  energy  can  be  expressed  as  a sum  of 
energies  of  each  degree  of  freedom,  the  distribution  function  for 
Maxwellian  equilbrium  may  be  expressed  in  factored  form, 


-fVXW,  ^(a'.^iv'K'CvO  -•  - CVj ) (2.i2) 


and  the  number  per  unit  area  of  diffusely  reflected  molecules  takes 
the  form 
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This  separation  of  variables  also  occurs  in  the  integration  for 
the  energy  flux  per  unit  area  due  to  diffusely  reflected  molecules. 


Making  the  further  assumption  that  the  energy  associated  with 
each  degreee  of  freedom  has  a quadratic  dependence  on  the  velocity 
involved,  a straight-forward  calculation  yields,  for  the  diffusely 
reflected  energy  flux  (Ref.  9,  p.  93) 
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For  the  non-Maxwellian  case,  one  cannot  factor  the  distribu- 
tion function  as  in  Eq.  2.12  and  the  corresponding  relations  for 
the  incident  energy  flux  presumably  would  become  very  complicated. 
As  a first  approximation,  however,  it  seems  reasonable  to  assume 
a similar  separation  of  translational  and  internal  energies  so  that 
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Combining  Eqs.  2.3,  2.9,  2.14,  and  2.15,  and, integrating 
over  the  cylinder  ^ ^ «, 

A = ™4k|  ff  [(ucN-v’+w'''  + jRT 

oc  ^ J J •*  J 
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3.0  RESULTS 
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Performing  the  integrations  indicated  in  Eq„  2.16  one  has 
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where  has  been  eliminated  by  the  relation 
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Jo  and  Li  are  the  modified  Bessel  functions  and  S is  the  molecular 
speed  ratio 
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(3.3) 


Eq.  3.1  is  valid  for  a monatomic  gas  ( j s O ) and,  under  the  assump- 
tions of  Eqs.  .2.14  and  2.15,  would  also  be  valid  for  a polyatomic  gas. 
The  results  of  Eq,  3.1  agree  with  those  of  Refs.  3 and  4 for  the  uni- 
form case.  It  will  be  observed  that  there  are  no  terms  corresponding 
to  a shear  stress  , nor  to  a heat  flux,  , in  the  transverse 

direction.  This  was  not  the  case  for  the  aerodynamic  force  calcu- 
lations of  Ref.  6.  Since  these  two  terms  are  the  principal  non-uni- 
formities in  a boundary  layer  flow,  it  follows  that  a heat  transfer 
type  of  free  molecule  probe  is  particularly  advantageous  for  boundary 
layer  surveys. 

If  > by  Ta*  , one  denotes  the  equilbrium  temperature  assumed 
by  the  cylinder  in  the  absenoe  of  heat  transfer  (adiabatic  wall  tem- 
perature), then  Eq.  3.1  may  be  solved  for  Taw/T  as  a function  , 
>p*K/p  J'Pi'j/f  an<*  In  the  resulting  expression  for  law  /T 

the  denominator  is  expanded  in  powers  of  ,-p^/p  and  « 

By  neglecting  squares  and  higher  powers,  one  obtains,  for  a monatomic 
gas. 


X 


aw 


= Hot 


+ 


4^u 

-f  th 


(3.4) 


T 


- 7 •» 


where i 


~ 'j'l  ^ “ 


Ff^  fo) 

G-o  (s) 

Gr.(S) 


'\ 


F.rs)  &«.,.(«) 

[S.fsf 

(Go  ("S  i} 


R « & 

&b(&) 


Res')  = (as1*  + 7s1  + 4)1.  (sVa)  +(2s‘'4-5sl)I1(s7iy 
1^(s)  = C2s’-+3/2)XofsVal  + (Zsl-i)x,(sVz) 

Ee,,(s\  = (sV2  f-3A)lo(sV,2)  -t-  (sV2+.  |)X,(S*/2) 

Fv(s)  = 2sl  I.  (s'/*}  + (i2sV£'Jr,(sVz) 

G.(s\  = (4s1  <-^I0(s»/2)  +•  4slI,(sVz) 

Qf»*(sW  I„(sV2^  - I,  (sVz) 

Sfw  («)  = X.  CWa)  + I,  ( SV2 ) 

&v  (s')  = - (‘tsl/5)Io('sl/a)  - (Hs2/£)X,  (&V2)  J 


The  functions  and  K^are  presented  graphi- 

cally in  Tig,  2. 

The  recovery  factor  is  defined  ly 
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where  "To  is  the  stagnation  temperature.  Using  the  relation 
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one  obtains 


a.  - 


Y 


Y-i 


or  using  Eq.  3.4 

1 


= 


Y-v 


I I Ta^j 
SM  r 


H.-l 


+ ^£aa_. 

S1 


Y> 

Eq.  3.9  may  be  written 
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These  "partial  recovery  factors*  are  presented  graphically,  for  a mona- 
tomic gas,  in  Fig.  3. 

One  defines  the  heat  transfer  coefficient  per  unit  area  as 
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Combining  Eqs.  3.1,  3*11  and  3*12 
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Eq.  3.13  may  be  written 
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These  “partial  Stanton  numbers*  are  presented  graphically, 
for  a monatomic  gas  in  Fig.  4. 

Inspection  of  Fig.  4»  reveals  that  the  partial  Stanton  numbers 
corresponding  to  the  non-unifonn  terms  are  all  small  compared  to  the 
Stanton  number,  $T©,  for  the  uniform  case8  Unless  a gas  flow  is  very 
non-uniform,  i.e.,  if  the  flow  itself  is  in  the  continuum  or  slip-flow 
region,  the  quantities  p^and  are  small  compared  to-p  . It  follows 
that  the  energy  balance  will  be  effected  only  slightly  by  the  contribu- 
tions arising  from  the  viscous  stress  terms. 
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